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ABSTRACT 

We present mapping results in several rotational transitions of HC3N, 
CeH, both cyclic and linear C3H2 and C3H, towards the cyanopolyyne peak 
of the filamentary dense cloud TMC-1 using the IRAM 30m and MPIfR 
100m telescopes. The spatial distribution of the cumulene carbon chain 
propadienylidene H2CCC (hereafter /-C3H2 ) is found to deviate significantly 
from the distributions of the cyclic isomer C-C3H2, HC3N, and CqB. which in 
turn look very similar. The cyclic over linear abundance ratio of C3H2 increases 
by a factor of 3 across the filament, with a value of 28 at the cyanopolyyne 
peak. This abundance ratio is an order of magnitude larger than the range 
(3 to 5) we observed in the diffuse interstellar medium. The cyclic over linear 
abundance ratio of C3H also varies by ~ 2.5 in TMC-1, reaching a maximum 
value (13) close to the cyanopolyyne peak. These behaviors might be related 
to competitive processes between ion-neutral and neutral-neutral reactions for 
cyclic and linear species. 

Subject headings: ISM: abundances — ISM: individual (TMC-1) — 
ISM: molecules — molecular processes 
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1. Introduction 

Among the molecules discovered in the interstellar medium, C3H2 and C3H are of 
peculiar interest for astrochemistry since both are observed in two isomeric forms: cyclic 
and linear. The ring molecule cyclopropenylidene (hereafter C-C3H2 ) was discovered in 
1985 (Thaddeus, Vrtilek, & Gottlieb) and has immediatly retained attention due to its 
ubiquity in the galaxy (e.g., Matthews and Irvine 1985). One of its linear counterparts, the 
propadienyhdene (Z-C3H2), was discovered in TMC-1 by Cernicharo et al. (1991). The two 
isomeric shapes of C3H, cyclopropynylidyne (C-C3H) and propynylidyne (/-C3H), were also 
detected in TMC-1 (Yamamoto et al. 1987; Thaddeus et al. 1985). 

This variety of isomers for the same species raises the question of their formation. In 
particular, are rings and chains formed from the same progenitors and involved in the same 
reactions networks? Furthermore these isomers can be used to probe interstellar chemistry 
models which include heavy interstellar molecules such as cumulene carbenes (Bettens and 
Herbst 1996, 1997; see also Millar, Leung, and Herbst 1987). Previous observations by 
Cernicharo et al. (1999) have shown that the cyclic over linear abundance ratio of C3H2 
(hereafter R2) in the diffuse medium along the line of sight toward the continuum sources 
W51E1/E2, W51D and W49 is one order of magnitude smaller than its value in TMC-1. 
In a recent study. Turner, Herbst, & Terzieva (2000) have compared measurements in 
three translucent clouds and two dense clouds (TMC-1 and L183) of the C3H cychc over 
linear abundance ratio (Ri) and R2 ratio - and column densities of 6 other hydrocarbons 
- to a modified version of the New Standard Model of chemistry (Lee et al. 1996). They 
found small variation of Ri from source to source, and a systematic higher R2 ratio with 
largest the values found in the two dark clouds. Also, it has been suggested that cyclic 
and linear isomers of a same compound could have different behavior in neutral-neutral 
and ion-neutral reactions (see, for example, Adams & Smith 1987; Maluendes, McLean, & 
Herbst 1993; Kaiser et al. 1997, 1999). The cyclic over hnear abundance ratio of molecules 
like C3H2 and C3H could then be used as a tool to investigate physical conditions in several 
media, from cold dark cloud to warmer and lower density media. 

In order to investigate these questions, we mapped in detail the region around the 
cyanopolyyne peak (hereafter CP) in TMC-1 in both /-C3H2 and C-C3H2, and observed at 
selected positions the two isomeric forms of C3H. CgH and HC3N have also been observed 
for comparison purpose. 
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2. Observations 

The C3H2 and CgH observations were made in 1990 and 1992 with the 100-meter 
telescope of the Max-Planck Institut fiir Radioastronomie at Effelsberg (Germany). We 
observed the lio — loi transitions of C-C3H2 at 18343.145 MHz, and simultaneously those of 
CeH^Hs/s J = 15/2 - 13/2 near 20794 MHz and /-C3H2 loi - Oqo at 20792.590 MHz. The 
antenna beamsize at 18.3 and 20.7 GHz are of 54" and 48", respectively. We used the 1024 
channel autocorrelator to achieve a resolution of 0.05 kms~^ at 18.3 GHz and 0.09 kms~^ 
at 20.7 GHz, and the data were taken by position switching the telescope. 

The HC3N and G3H observations were done in 1995 and 1999 with the IRAM 30-meter 
telescope at Pico Veleta (Spain). We observed the I-C3H ^H3/2 J = 9/2 — 7/2 transition 
near 98 GHz and at the same time the transition of C-C3H 2i2 — In near 91.5 GHz. The 
observations were made in frequency-switching mode. The autocorrelator was used as a 
spectral instrument with a velocity resolutions of 0.12 kms"^. The half-power beamwidth 
(HPBW) and main-beam efficiency are 26" and 0.75 for /-C3H and 27" and 0.78 for C-C3H. 
System temperatures were in the range 90-130K. For HC3N, observed at 90979.023 MHz 
(J = 10 — 9), the HPBW is 26". Pointing and calibration were monitored by regularly 
observing planets and quasars for both telescopes. Sample spectra are shown in Fig.l and 
line parameters for selected positions are given in Table 1. 

3. Results 

3.1. C3H2 isomers and CeH 

The mapping results of this study are shown in Fig. 2. The distributions towards 
the CP in TMC-1 of C-C3H2 and 1-Csli2 are compared to that of CeH in Fig. 2a and b, 

respectively. Fig. 2c shows the distributions of CgH and HC3N. The spatial distributions of 
/-C3H2 and C-C3H2 (Fig. 2d) are clearly different, whereas the emissions of C-C3H2, CgH 
and HC3N have comparable distributions. The emission of /-C3H2 is shifted toward the 
west by approximatively 40" with respect to C-C3H2 or CeH - note that, since CeH and 
/-C3H2 were observed in the same bandwidth, this shift cannot be due to pointing errors 
during the observations. 

In order to derive relative abundances, we first computed relations between the 
column-density and the observed line intensity valid for uniform physical conditions, i.e. 
n(H2) = 3x10^ cm~3 ^ ^qK in TMC-1 (Cernicharo & Cuclin 1987). From a 

statistical equilibrium calculation using coUisional excitation rates for C-C3H2 from Avery 
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& Green (1989), and assuming an ortho/para ratio of 3, we find: A^(c-C3H2)[cm~^] = 
2.2x10^^ /Tmfefif [K-km-s~^] . This linear relationship indicates that the lio — loi transition 
of C-C3H2 is optically thin in TMC-1; it is in good agreement with the work of Cox et 
al. (1989). Madden et al. (1986), from observations of the isotopic c-C2^^CH2 in TMC-1, 
derive an optical depth ranging from 5.6 to 6.8 for the 2i2 — loi transition of C-C3H2 at 85 
GHz. This result is in agreement with ours since the 2i2 — Iqi line has a larger opacity 
than the lio — loi transition for dark cloud physical conditions. Indeed, with our code, we 
are able to reproduce the results of Madden et al. at 85 GHz, the lio — loi transition still 
being thin or marginally saturated (r < 1.5). For /-C3H2 and CeH there are no coUisional 
cross sections available. We have estimated cross sections from those of H2CO (calculated 
by Green et al. 1978) and of HC3N (Green & Chapman 1978) - sec Ccrnicharo et al. 
(1999). This gives: iV(/-C3H2)[cm-2] = 2.2x10^3 J T^,^dv[K-km-s-^] and A^(C6H)[cm-2] = 
5.2x10"'^^ f Tmbdv[K-kTa-s^^] . Note that for CeH, JTmbdv is the total integrated intensity 
summed over the four components of the hypcrfinc structure. These relationships are in 
good agreement with previous works (Ccrnicharo et al. 1991; Bell ct al. 1999). 

Comparing two positions separated by 40" (which corresponds to 0.02 pc for the 
adopted distance of 100 pc to TMC-1 - Ccrnicharo & Guelin 1987), we derive the following 
column density: at the CP (0,0), 7V(C6H) = 8.3x10^2 ^^^-2^ 7V(c-C3H2) = 5.8x10^3 ^^^-2^ 
iV(/-C3H2) = 2.1x10^2 (R2=28±6); and at the edge of the TMC-1 filament (-40,0), 

NiCeii) = 4.7x10^2 cm'^, A^(c-C3H2) = 2.8x10^^ cm'^, A^(;-C3H2) = 3.2x10^^ cm^^ 
(R2=10±3). Between these two positions, the cyclic over linear abundance ratio changes 
by a factor of ~ 3. This variation cannot be an artefact caused by calibration errors. 
Indeed, while the calculated CeH over /-C3H2 column density ratio varies from 1.5 to 4, the 
C-C3H2 over CgH column density ratio remains constant {c^ 6.5). Since the former ratio is 
unaffected by calibration and pointing errors (lines arc observed in the same bandwidth), 
the observed variation of R2 must be real. This result suggests that chemical gradients are 
present in TMC-1 on scales smaller than 0.02 pc. 

To investigate further the spatial variations of molecular abundances across the TMC-1 
filament, we averaged the spectra along the six different rows shown in Fig. 2. Results are 
summarized in Fig. 3: R2 ranges from 12 to 37 and clearly increases from the west to the 
east of TMC-1. Note that a constant density of n(H2) = 3x10"^ cm^^ has been adopted. If 
we use instead H2 densities derived by Pratap et al. (1997) from an analysis of the HC3N 
transitions while keeping the same kinetic temperature (10 K), the values of R2 are then 
lowered by 20 to 40 % but the west to east variation of R2 remains. The rise of the cyclic 
over linear abundance ratio of the C3H2 isomers at a spatial scale of 2 arcmin (0.06 pc) is 
therefore a firm result. 
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3.2. C3H isomers 

As C3H is thought to be formed by the same reaction as C3H2, i.e. by the dissociative 
recombination of C3H3" (Adams & Smith 1987), it is interesting to compare N^c-C-^H) / N{1- 
C3H) (hereafter Ri) with R2. What is the value of Ri and does it vary across the filament 
similarly to R2? In order to answer these questions we observed C-C3H and /-C3H at 
13 positions along two strips close to the CP. One strip crosses the filament from offset 
positions (60" ,60") to (-60" ,-60"), the other from (10,-70") to (-70", 10"). We also observed 
the C-C3H2 and the /-C3H2 peaks. Where possible, we analyzed the data with the HFS 
method of CLASS (a software developed by the GILDAS working group). This method 
provides the total optical depth, the average linewidth and the brightness temperature of 
a line with hyperfine structure. A reliable estimate of the excitation temperature could 
be obtained for several points. We found that 3K < T^^ < 3.8K for C-C3H (in excellent 
agreement with Mangum & Wootten 1990), and 4.9K < T^o, < 6.7K for I-C3R. In the 
following, we adopt Tex== 3.5K for C-C3H and Tex= 5.5K for the hnear isomer. 

Calculations of the column densities were done using the classical formula with n 
= 2.4 D (Yamamoto et al. 1987) for C-C3H and = 3.1 D (Green 1980) for /-C3H. For 
each observed point, we computed the average of the total column densities by using each 
hyperfine transition weighted by the inverse of the square of the error. The resulting column 
densities and the values for Ri are given in Table 2. Comparing with the values of R2 in 
Fig. 3, we see that: 1) With a value of 12 at the CP, Ri is smaller than R2; 2) As R2 does, 
Ri shows variations across the filament (by a factor of 2.5), mainly due to variations of 
Nic-CsR). 

4. Discussion 

The observations described above underscore the differences in the distribution and 
relative abundances of the cyclic and linear forms of C3H2 and C3H in the molecular 
filament TMC-1. In the following, we study the origin of this behavior and suggest that it 
is driven by competitive processes between ion-neutral and neutral-neutral reactions. 

4.1. Steady-state calculation 

We assume here, as it has been proposed (see, for example, Adams & Smith 1987; 
Maluendes, McLean & Herbst 1993), that cyclic C3H2 and C3H are both formed in the 
dissociative recombination of C-C3H3 , while linear C3H2 and C3H result from H2C3H+. We 
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also assume that the C3H+ + H2 association reaction is the dominant formation mechanism 
of C3H3 , which is produced in equal amounts in cychc and hnear forms. 

At steady state, the chemical kinetic equation for the abundance of C-C3H2 - that is, 
x(c-C3H2) - is: 

dx(c-C3H2)/dt = 2b^krx(c-C3H+)Xe - 2Kfx(c-C3H2) = 

where "b" is for "branching ratio", is for "recombination rate" [cm^-s~^], K'^ is for "mean 
destruction rate" [s~^] (K'^ = J2ikf ■ Xi , where Xj is the second reactant). We use indices 
(l,c) and exponants (2,1) to distinguish between C3H2 and C3H in linear and cyclic forms 
(for example, ^bc means "branching ratio for cyclic C3H2"). Similar equations can be 
written for /-C3H2, C-C3H and /-C3H. Using the observed column densities for each species 
at the border of the filament (i.e. row number 2), where the electron density should be the 
highest in our data set, and writing k^ = k[ and X = x(c-C3H3 )/x(l-C3H3"), calculations 
lead to: 

x(c-C3H2)/x(c-C3H) = (2b,/ibe)(iK^/2K^) ~ 5 ; 
x(/-C3H2)/x(/-C3H) = {^h/%)CKfpKf) ^ 2 ; 
R2 = (%A,)(2Kf/2K^)X~16; 
Ri = (VbO('Kf/iK^)X~6. 

A simple and coherent set of solution for this system is ^b; ~ ^b; ~ ^bc — ^bc ; ^Kf ~ 2 ^Kf 
; ^K^ ~ 5 ; ^Kf ~ (1/6) ^Kf X ; ^K^ ~ (1/16) ^Kf X . 

Under the hypothesis of steady state and from the observed abundance ratio at the 
border of the TMC-1 filament, we find that C-C3H2 is destroyed approximatively 5 times 
slower than C-C3H and that the /-C3H mean destruction rate is twice the /-C3H2 one. 

To our knowledge, no additional formation process of C-C3H2 in dark clouds can be 
invoked to explain an increase of its abundance with respect to /-C3H2. Variations of R2 
(and Ri) in TMC-1 could result from variations of the destruction rates of C-C3H2 and 
/-C3H2. Such an explanation has been proposed by Maluendes, McLean, & Herbst (1993): 
while C-C3H2 is inert with respect to most neutral-neutral reactions, /-C3H2, as /-C3H and 
C-C3H, arc easily destroyed by these reactions. Hence, because neutral-neutral reactions 
proceed fast when a radical reacts with abundant atoms, it is possible in a dark cloud like 
TMC-1 that R2 is one order of magnitude larger than in a more diffuse medium where 
the proportion of reactive ions (increasing with respect to reactive atoms) is sufficient 
to destroy C-C3H2 and ^-C3H2 at the same rate. The same explanation holds for the 
progenitor ions c-CsH^ and I-C3HJ (see Cernicharo et al. 1999): while I-C3HJ - assumed 
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to be the progenitor of /-C3H2 - can be efficiently removed through ion-neutral reactions 
lowering the amount available to produce /-C3H2, c-CsH^ is mainly affected by dissociative 
recombination to produce C-C3H2. R2 thus depends on the ion-neutral reactions of l-CsHg" 
that do not affect the cyclic ion: in dark clouds, where the abundance of reactive molecules 
can be hundred times larger than that of electrons, I-Caffjj" is removed faster by ion-neutral 
reactions than by the dissociative recombination lowering the production rate of /-C3H2, 
and increasing R2. In a more diffuse medium, where the ion-neutral reactions are much 
less efficient, I-C3H3 and C-C3H3" are removed by the same process and R2 is then closer 
to ^Kf/^Kf (whose value also decreases with decreasing density as suggested above). This 
explanation is strengthened by the lower R2 values found in the diffuse medium from 
absorption measurements by Cernicharo et al. (1999) - 3 to 5 versus 10 to 40 in TMC-1 - 
and by the observed Ri variations across the filament. We test this hypothesis in the next 
section. 



4.2. Chemical modelling 

In order to study the impact of physical conditions on the R2 ratio, we have run 25 
models of gas phase chemistry with different densities and visual extinctions using a time 
dependent chemistry code solving the system of stiff ordinary differential equations with 
the Gear method. Although we have used a time dependent chemistry code - samples of 
the evolution chemistry arc shown in Fig.4 - we discuss here the results at steady state. 
Indeed, we would like to study the basic processes of the isomeric differentiation, and to 
identify the dominant reactions leading to the observed large variations of R2 in diffuse and 
dense clouds. We aim first at a qualitative description of the mechanism. A quantitative 
description would require a more accurate knowledge of the reaction rates (including 
branching ratios) important for C-C3H2 and /-C3H2. Fig.4 shows that C-C3H2//-C3H2 and 
c-C3H3'/l-C3H^ abundance ratios reach values close to the steady state value, just after 
they reach their maximum abundance. As a first step, it is thus reasonable to consider 
only the processes leading to variation in isomeric abundances at steady state. We have 
used, as chemical network, the UM1ST95 database (Millar, Farquhar, & Willacy 1997) in a 
pure gas-phase scheme. As recommanded by the authors, the species C2H5, C2H6 and their 
ions have been excluded from the network. We have also excluded the species including 
the following atoms: P, Si, CI, Na, Mg. On the other hand, we have included the grain 
surface formation of H2 with a rate coefficient 1.5x10"^^ n^(H) cm~^ s~^. Moreover, we 
have updated the dissociative recombination rate of C-C3H3 at 300 K to 7.0±2.0 x 10~^ 
cm^ g-i (Abouelaziz et al.l993) and assumed the same value for I-C3H3". Parameters of the 
models and initial elemental abundances are given in Table 3 and 4, respectively. 
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In Fig. 4, we see that R2, and also X, increase with increasing visual extinction. This 
general trend of our models agrees with observations since R2 is lower in the diffuse medium 
than in the TMC-1 dark cloud. The models also show that the increase is faster in higher 
density media. We have plotted in Fig. 5 the electronic density with respect to R2 for each 
model. The distribution of the points exhibit clearly a correlation between the ionization 
fraction and the cyclic over linear ratio which can be separated in two regimes: "low" 
{l<Av<2) and "high" (5<Ay<10), Ay = 3 mag being an intermediate case. In the low 
extinction regime, R2 strongly depends on the electronic abundance. This is an interesting 
result because it opens the possibility to use R2 as a tool to probe the electronic abundance 
in low extinction regions where H^^CO+ and DCO"*" cannot be detected. Why is the cyclic 
over linear ratio sensitive to the ionization fractionation? In the UMIST95 chemistry, C"*" 
can destroy both /-C3H2 and C-C3H2, whereas C reacts only with the linear isomer. R2 
is then sensitive to the electronic abundance via the C"*'/C ratio as shown in Fig. 6. Note 
that the range 3 to 5 that Cernicharo ct al. (1999) have observed in the diffuse medium 
corresponds, as expected, to high electronic abundances and to moderate visual extinctions. 
In the high extinction regime, R2 is independent of the C+/C ratio, while it is sensitive to 
the (H3 +C+)/0 ratio (Fig.7). Indeed, in the UMIST95 database the atomic oxygen reacts 
only with Z-C3H2 (O + /-C3H2 — > C2H2 + CO). This reaction is neghgible with respect to 
other destruction reactions (in particular the reactions with C"*" which affects both isomers) 
until all the carbon is locked into CO. It is thus unefficient in low visual extinction media 
where the abundances of C and C"*" remain high. On the other hand, it becomes the 
principal destruction channel of /-C3H2 at high visual extinctions: for Ay = 2 mag and a 
density of 10^ cm~^ the O + /-C3H2 reaction is roughly 15% of the /-C3H2 destruction rate 
while it represents more than 80% of it at Ay — 10 mag. The variations of the cyclic over 
linear abundance ratio of C3H2 can thus be understood as a consequence of the competition 
between neutral-neutral and ion-neutral reactions in the interstellar medium. 

The determination of electronic abundances from measurements of R2 using Fig. 5 may 
not be straightforward. Indeed, taking into account the large uncertainties in the reaction 
rates and branching ratios, it is difficult to infer quantitative results from this study on 
the C3H2 isomers. We can say however that low R2 values are specific of high electronic 
abundance while high R2 values indicate low electronic abundance. In a more general 
way, isomeric ratio can be used to probe physical conditions in different media from the 
diffuse gas to dark clouds. In the future, cyclic and linear species other than the isomers of 
C3H2 shall be included in chemical models. C3H is a good candidate. Indeed, Kaiser et al. 
(1997, 1999) have shown that the neutral-neutral reaction C + C2H2 — > C3H + H produces 
linear and cyclic C3H, with an increase of the linear over cyclic abundance ratio with rising 
collision energies. It would be interesting to make measurements of Ri in absorption in 
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the diffuse medium in order to compare this neutral-neutral formation path with the C3H3" 
dissociative recombination. 

This study has been partly begun. Indeed, in a recent paper, Turner, Herbst, & 
Terzieva (2000) investigate the Ri and R2 ratios in three translucent clouds and in the 
dark clouds TMC-1 and L183. They compare observations to the predictions of a modified 
version of the New Standard Model of chemistry (Lee et al.l996) which includes cyclic and 
linear species for all precursor to C3H and C3H2. One of their conclusion is that the clear 
observed anticorrelation of R2 with extinction suggests that the ratios are strongly affected 
by formation or destruction rates, and not just by branching ratios among the relevant 
chemical reactions. 



5. Conclusion 

In this paper, we have presented new results about the cyclic over linear abundance 
ratios of C3II2 and C3H toward the CP of TMC-1. These ratios exhibit variations which 
are probably due to variations in the electronic abundance across the filament. Indeed, the 
chemical modelling of the reactions that affect C-C3II2 and /-C3II2 shows that their ratio is 
sensitive mainly to the behavior of the C+/C (at low visual extinction) and (H3 +C+)/0 
(at high visual extinction) abundance ratios, which are both related to the electronic 
abundance. As C3H2 is observed in its cyclic and linear form in different physical conditions, 
it could be used as a tool to probe the electronic abundance in the ISM. However, to infer 
accurate values of the fractional ionization from cyclic and linear isomers abundances, the 
branching ratios and kinetic rates used in chemical networks must be more accurately 
known. In addition to laboratory studies, a systematic measurement of Ri and R2 in 
different astrophysical media should be useful for this work. 
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Table 1: Line Parameters for CgH, Z-C3H2 and C-C3H2 at Selected Positions in TMC-1 
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5.89±0.02 


0.20±0.02 


















2.46±0.59 


6.12±0.02 


0.14±0.03 


-40,0 


0.31±0.03 


6.74±0.01 


0.2B±0.00 


0.31±0.07 


6.78±0.03 


0.44±0.08 


-50,0 


2.20±0.76 


6.48±0.06 


0.20±0.00 




O.lliO.ll 


5.99±0.02 


0.09±0.20 










3.99±0.76 


5.70±0.03 


0.20±0.00 


0,40 


0.29±0.04 


5.82±0.03 


0.21±0.06 


0.28±0.06 


5.85±0.02 


0.18±0.04 


0,50 


3.96±0.71 


5.68±0.02 


0.17±0.04 




0.25±0.04 


6.07±0.04 


0.21±0.09 


0.26±0.06 


6.15±0.02 


0.14±0.04 




4.58±0.71 


5.96±0.02 


0.28±0.05 


40,-40 


0.48±0.03 


5.78±0.01 


0.25±0.00 


0.28±0.07 


5.80±0.03 


0.21±0.0fi 


50,-50 


2.83±0.78 


5.61±0.03 


0.19±0.05 




0.22±0.03 


6.03±0.03 


0.27±0.05 


0.23 = 0.07 


6.08±0.03 


0.13±0.0(i 




4.88±0.78 


6.88±0.03 


0.34±0.06 


Notes - 


Positions are 


offsets relative 


to the Cyanop 


olyync Peak : 


a(1950) = 


04''38'"38.6=, 


i5(1950) = 


25°35'45.0" 







TMC-1 exhibits two main components at vlsr — 5.3 and 6.1 km s ^ (Sume, Downes, &c Wilson 1975) 
When possible, we have fitted several components along a line of sight 



-14- 



Table 2: Derived Column Densities 



for C-C3H and /-C3H in TMC-1 



Position 




iV(/-C3H) 


Nc/M 


( " , " ) 


10^2 cm-2 


10^^ cm-2 


Ri) 


-40,-40 


3.7±0.2 


<4 


>9 


-20,-20 


7.7±0.3 


8.0±0.2 


10±0.5 


0,0 


10.3±0.3 


8.4±0.1 


12±0.5 


20,20 


10.0±0.3 


7.6±0.4 


13±1.2 


40,40 


5.8±0.3 


6.7±0.5 


9±1.1 


60,60 


3.6±0.2 


6.1±0.7 


6±1.0 


10,-70 


5.2±0.2 


7.0±0.4 


7±0.6 


-10,-50 


5.3±0.2 


7.8±0.2 


7±0.4 


-30,-30 


5.1±0.2 


8.0±0.1 


6±0.4 


-50,-10 


4.7±0.2 


7.1±0.5 


7±0.8 


-70,10 


3.7±0.2 


7.3±0.3 


5±0.5 


30,-60^ 


7.8±0.3 


8.0±0.2 


10±0.6 


-30,0^ 


10.0±0.4 


8.2±0.2 


12±0.7 


Notes - Positions are 


offsets relative to the CP : 


a(1950) ■- 


= 04'^38"*38.6^ (5(1950) = 


25°35'45.0" 


^ C-C3H2 


peak 






2 Z-C3H2 


peak 
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Table 3: Parameters of the Chemical Models 



Parameter 






Value 


n(H2) (103 


cm 3) 


1, 


3, 10, 30, 100 


T(K) 






10 


Ay (mag) 




1 


, 2 , 3 , 5 , 10 


C (s-^) 






1.3 X 10-^^ 


iLi (erg cm" 






2 X 10-15 



C is the cosmic-ray ionization rate 

Ui is the density of ionizing radiation at Ay = 
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Table 4: Initial Elemental Abundances 
with respect to Hydrogen 



He 1.4 X 10-1 

C 7.3 X 10-5 

N 2.14 X 10-5 

O 1.76 X 10-^ 

S 2.0 X 10-^ 

Fe 3.0 X 10-9 
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Figure Captions 

Fig. 1. — Spectra of CeH, cyclic and linear C3H2, and cyclic and linear C3H towards the 
cyanopolyyne peak of TMC-1 {vlsr = 5.8 kms"^). Note that the CqH ^Hg/a J=15/2-13/2 
quadruplet and the ^-C3H2 loi-Ooo transition are observed in the same bandwidth. The 
J=3/2-l/2 F=2-l and F=l-0 transitions of C-C3H 2i2-lii are not shown here. 

Fig. 2. — Contour levels of integrated intensity in TMC-1 for: Panel (a): CeH ^113/2 
J=15/2-13/2 and C-C3H2 lio-loi; Panel (b): CeH ^U3/2 J=15/2-13/2 and /-C3H2 loi-Ooo 
(the crosses indicate the measured positions); Panel (c): CeH ^113/2 J=15/2-13/2 and HC3N 
10-9; Panel (d): C-C3H2 lio — loi and /-C3H2 loi — Oqo (numbers indicate the values of R2, 
i.e. the ratio A^(c-C3H2)/A^(/-C3H2)). All contours go from 50 to 100 % of the peak value 
by steps of 10 %. The peak values are 0.16 K-kms^^, 2.6 K-kms^^, 0.15 K-kms^^, and 1.75 
K-kms^^ for CgH, C-C3H2, /-C3H2, and HC3N respectively. The coordinates of the (0,0) 
position arc: a = 04^38™'38.6'^, 6 = 25°35'45.0" (1950.0) corresponding to the cyanopolyyne 
peak. The boxes to the right compare the CeH and /-C3H2 spectra averaged along the six 
directions indicated by the arrows. Note the significant emission of Z-C3H2 toward the west 
side of the TMC-1 filament where the emission of CeH decreases rapidly. 

Fig. 3.— Variations of R2 (i.e. the ratio iV(c-C3H2)/A^(/-C3H2)) across the TMC-1 filament 
using: Lower panel: a constant density of 3x10"^ cm~^; Upper panel: densities derived 
by Pratap et al. (1997). In order to increase the signal to noise ratio, the spectra have been 
averaged along six rows (separated by 28") parallel to the TMC-1 filament. Lower limits in 
row number 6 are derived from 3a upper limits on the /-C3H2 emission. The cyanopolyyne 
peak (CP) is located within row number 4. 

Fig. 4. — Chemical evolution for three different visual extinctions at n(H2) = 10^ cm~^'^ and 
T = 10 K. Upper boxes show abundances relative to H for C, /-C3H2, C-C3H2 and e~. Lower 
boxes show the resulting cyclic over linear abundance ratio for C3H2 and C3H3 isomers. 

Fig. 5. — The predicted electronic abundance vs. R2 at steady state. Two regimes can be 
distinguished according as the visual extinction is low or high. 

Fig. 6. — The predicted dependence of R2 with the C+/C ratio at steady state. A clear 
relationship is established at low visual extinction between R2 and the ionization fraction of 
atomic carbon. Note that at high visual extinctions R2 depends only weakly on C+/C. 

Fig. 7. — Model prediction for the relation of R2 vs. the (H3 -|-C+)/0 ratio at steady state. 
Note the good correlation between R2 and (H^-|-C''")/0 at high visual extinctions. 
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